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Abstract 



S-h ' We present a model for quark matter with a density dependent quark- 

quark (confining) potential, which allows to describe a deconfinement 
phase transition as the system evolves from a low density assembly of 
bound structures to a high density free Fermi gas of quarks. A proper 
account of the many-body correlations induced by the medium is cru- 
cial in order to disentangle this behaviour, which does not uniquely 
stem from the the naive density dependence of the interaction. We 
focus here on the effects of finite (non-zero) temperatures on the dy- 
namical behaviour of the system. In particular we investigate the 
ground state energy per particle and the pair correlation function, 
from which one can extract the relevant information about the quarks 
being confined or not; the temperature dependence of the transition 
density is also derived. 
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1 Introduction 



"QCD inspired" models for quark/nuclear matter have been proposed since 
many years, mainly within the framework of non-relativistic constituent 
quark models |], @, |3], |fl; the latter have proven remarkably useful for single- 
hadron spectroscopy, in spite of the lack of Lorentz invariance and chiral sym- 
metry. Relativistic models, based on effective chirally invariant Lagrangians 
of NJL type or related to the so-called colour dielectric model, have also 
provided important tools for the investigation of quark matter |5|, |5], [7|, ^ [J; 
however their many-body description mainly relies on the mean field ap- 
proach, which does not seem the most appropriate one to deal with the 
problem of a quark-hadron phase transition. 

The present work is based on a non-relativistic model for a one-dimensional 
many quark system, whose key feature is a density dependent interaction be- 
tween quarks [§(]: 

V(x) = \x 2 e~ cp W , (1.1) 

x = x\ — %2 being the relative interquark distance, p = N/L the (uniform) 
density of the system with N fermions in a length L, and c a constant pa- 
rameter. The potential Ql.lj ) resembles a strong confining force in the limit 
of very low densities, where one expects quarks to be bound into hadrons, 
while it becomes negligible at large densities, where the quarks behave as a 
free Fermi gas (which should mimic the quark-gluon plasma phase). 

Internal degrees of freedom are neglected, but antisymmetry for the global 
wavefunction of the system is required, as it is appropriate for an assem- 
bly of fermions. We notice that while the many-body interaction potential 
of 0, can only be dealt with variational Monte Carlo techniques, the 
present two-body potential allows to describe the system within the cus- 
tomary many-body schemes, which have been widely tested against nuclear 
matter properties. 

The philosophy underlying this model is to provide a tool for a phe- 
nomenological but microscopic description of the phase transition from a 
plasma of (weakly interacting) quarks and gluons to an assembly of colour- 
less quark clusters, namely hadrons. This transition should have occurred in 
the early stages of the Universe, at rather large temperature and density of 
the primordial gas of elementary constituents; as the system expanded and 
cooled down, hadronization took place, giving rise to the ordinary matter of 
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baryons and mesons, eventually bound into nuclear systems. 

Relativistic heavy ion collisions are presently investigated in an attempt 
of reproducing the extreme density and temperature conditions at which the 
quark-gluon plasma phase can occur. One can thus envisage that tempera- 
ture plays a crucial role in the description of the dynamical evolution of the 
system. Indeed, according to lattice QCD calculations, confining forces be- 
come weaker with increasing temperature, as it is suggested by the following 



temperature dependence of the string tension |10 
cr(T) = cr 



d ( T c-T), (1-2) 

where <7o = cr(T = 0) and the theta function implies that confinement no 
longer survives above the critical temperature Tq- Obviously this can only 
be taken as a qualitative indication of how the temperature affects the quark- 
quark interaction potential; moreover various uncertainties concern the pre- 
cise determination of the critical temperature (as well as the order of the 
phase transition which shows up within this framework). Nevertheless we 
have assumed Q1.2 ) as a sensible phenomenological Ansatz for the tempera- 



ture dependence of the quark-quark interaction strength and employed it in 
our model potential (1.1) in order to describe the dynamical evolution of the 



system at finite temperature. 

In Section 2 we shall shortly review the many-body formalism employed 
to deal with the system at zero temperature, introducing the relevant quan- 
tities which might signal, as a function of the density, the occurrence of 
the above mentioned phase transition. Numerical results are reported at 
T = in order to better appreciate the modifications introduced by finite 
values of the temperature. The latter are explored in Section 3, where a few 
schematic details of finite temperature field theory precede the description of 
the many-quark dynamics at finite temperature. The results for the ground 
state energy per particle and for the pair correlation function are presented at 
various (increasing) temperatures within the framework of the independent 
pair approximation. Finally Section 4 summarizes the virtues and limitations 
of the present approach. 
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2 Results in the T = limit 



Before introducing the finite temperature description, we shall report here 
an outline of the results obtained at zero temperature. We have mainly 
investigated the ground state energy (per particle) and the pair correlation 
function at different values of the density. 

The medium induced correlations have been taken into account by solving 
the Bethe-Goldstone equation for the wavefunction of an interacting pair, 
within the independent pair approximation. For this purpose we expand the 
(antisymmetrized) two-particle wavefunction as follows: 

e iKX ( \ 

= N—j=- < c n0 sin(fcx) + c nj sm(kjx) > , (2.1) 

VL [ k l3 ,k 2] >F J 

where X — {x\ + x 2 )/2, x = X\ — x 2 , kj = (kij — k 2 j)/2 and the total 
momentum K = kij + k 2 j is conserved in all terms in the expansion. The 
restriction on k\j and k 2 j to be above the Fermi momentum kp accounts for 
the Pauli blocking of the states available to the considered pair, which is due 
to the presence of the medium. In the present one-dimensional calculation 
it can be rewritten as: 

kj > k F + \K\/2. (2.2) 

The trial wavefunction ( |2.f|) is then inserted into the Schrodinger equation 
for the relative motion with the potential ( |1.1| ): 

i) n k{x) = E nk i/j nk (x), (2.3) 

which, by exploiting the orthogonality of the unperturbed wavefunctions in 
( |2.1| ), reduces to a system of linear equations for the coefficients c n j of the 
Bethe-Goldstone solution we are looking for. One should also keep in mind 
that we are interested in bound (or quasi-bound) states of the interacting 
pair, at least for those (low) densities where the confining potential ( |1 . 1| ) is 
strong enough to actually provide confinement. In order to disentangle the 
few discrete eigenvalues of ( |2.3j ) out of its continuum spectrum, we impose 
on ip n k(x) the boundary condition 

ipnk(R) = 0, (2.4) 



d 2 



dx 2 



+ V{x) 
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R being an arbitrary length, of the order of (or larger than) the peak position 
of the potential barrier developed by (|1.1|) at intermediate distances. Eq. 
( |2.4p is easily satisfied by choosing the wavenumbers in the expansion ( |2.1| ) 
to be integer multiples of ir/R. If ip n k corresponds to a truly bound state, 
then its vanishing at \x\ = R will not be accidental and we expect such a 
solution (and the corresponding energy eigenvalues) to be fairly stable with 
respect to broad variations of R. 



At zero temperature and very low densities (p < 0.1) we have found |TT 
that the Bethe-Goldstone wavefunctions are close to the bound states one 
can obtain in the pure Schrodinger equation (namely in the absence of Pauli 
blocking), although they develop some small components outside the "con- 
finement" region. The importance of these components rapidly increases with 
the density, thus showing that the effect of the medium on the interacting 
pair loosens its binding and produces a gradual transition (at intermediate 
densities) to a phase of non-interacting particles (Fermi gas). 

By summing over the energy eigenvalues of particles in occupied states, 
namely those with relative momentum ki < kp — \K\/2, one can evaluate the 
ground state energy of the system: its evolution as a function of the density 
is shown in Fig. 1, where it is compared with the purely kinetic energy of a 
Fermi gas. The two coincide for large densities (say p > 1.0) but in the low 
density range the effect of correlations allows to disentangle the existence 
of dynamically bound pairs, whose energy is fairly independent upon the 
density of the system. 

Another quantity which provides a clear signature for the existence (if 
any) of bound clusters in the system is the so-called pair correlation function 
(r = \x\) 

p 2 

which is obtained as the ground state expectation value of the two-body 
density operator: 

p 2 (x 1 -X 2 ) = J2 S ( X i ~ X M X j - X 2)- ( 2 - 6 ) 

Within the independent pair approximation g(r) turns out to be: 

g B c(r) = 4r E IVW*)I 2 [k = (k 1 -k 2 )/2] (2.7) 
p L lfcll,lfc 2 l<fc F 
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Figure 1: The ground state energy per particle at T = calculated in the 
independent pair approximation (black dots) is shown as a func- 
tion of density, together with the one of a free Fermi gas (dashed 
line) . 



and can thus be evaluated using the ground state Bethe-Goldstone wave- 
functions (|2.1j). In Fig. 2 we display a few g BG (r) at T = and for different 
densities: one can appreciate the evolution of the pair correlation function 
from the low-density regime, where confinement is evident up to p ~ 0.15 
(lower densities, not displayed here, show even higher peaks at small dis- 
tances), to the transition density (p ~ 0.19), where the Fermi gas component 
starts developing at large distances, up to densities of the order of 0.5 and 
higher, where the pair correlation function practically coincides with the one 
of an uncorrelated Fermi gas. The correlation function obtained within the 
present model compares fairly well with the one of ref. jj, thus showing the 
equivalence of the two approaches for what concerns the physical properties 
of the system. 
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Figure 2: The pair correlation function g(r) derived from the Bethe- 
Goldstone wave functions at T = is displayed as a function 
of the relative interquark distance: the various curves are labelled 
by the value of the density. The free Fermi gas correlation func- 
tion is shown only for p = 0.5 (dotted line), where it practically 
coincides with the one of the correlated system. 

It is worth noticing that, in the presence of strong short range interactions, 
the importance of the medium grows with the density (which is proportional 
to &f); in the present model, however, the interparticle potential exponen- 
tially decreases with the density and thus the largest modifications on the 
relative motion of a pair of particles only occur within an intermediate, lim- 
ited range of densities. This is an interesting interplay between Pauli and 
dynamical correlations. 
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3 Finite temperatures 



In order to deal consistently with the hadron-quark plasma phase transition 
one should take into account the behaviour of the system with the temper- 
ature, thus extending the microscopic description within the formalism of 
finite temperature field theory. In this case it is convenient to statistically 
treat the system within the grand canonical ensemble, by defining 



K = H-fj,N, 



(3.1) 



N being the number operator and fi the chemical potential. 

The expectation value of any operator will then be evaluated by imple- 
menting the ensemble average 



< O >= Tr(p G 0), 

Pg being the statistical density matrix (f3 = 1/ksT) 

1 



Pg 



-(3K 



'G 



and Zq the grand partition function 



(3.2) 



(3.3) 



(3.4) 



In the above the trace (Tr) implies a sum over a complete set of states in the 
Hilbert space with any number of particles. 

Moreover one can introduce the so-called temperature Green's functions: 
for example the single particle propagator is defined as 



£ a/ j(xr, xV) = -Tr {p G T T tjj Ka (xT)ijj ] Kl3 (x.'T') } 



(3.5) 



where ^(xr) is the imaginary-time field operator in the (modified) Heisen- 
berg picture [JEJ and T T the corresponding time ordering operator. 

For a non-interacting system of fermions the temperature Green's func- 
tion reads then (in momentum space): 



g° a Jk,cu;T)=8, 



a/3 



n° k (T) 



+ 



10 



0° -p)+iv u - {el - fi) - %Tf 



(3.6) 
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where the occupation probability 

n° k (T) 



1 



1 4. 6 /3(^-m) 



(3.7) 



reduces to the sharp theta function momentum distribution of the Fermi gas 
in the limit of zero temperature: 



<{T) 



{el -el) 



(3.8) 



These ingredients can be used to evaluate, at T ^ 0, the polarization 
(particle-hole) propagator within the customary perturbation theory and, 
from it, to extract the two-body correlation function according to the pro- 
cedure described in ref. [|TTJ. 

It should be pointed out that, in principle, the chemical potential has to 
be self-consistently determined in order to satisfy the relation: 



N 



= P 



de 



y/2,71 Jo y/e[l + exp{(3(e - //)}] 



(3.9) 



which, in the limit of very small temperatures leads, for a free Fermi gas, to 
the expansion: 



LL~E F 




(3.10) 



The latter, however, becomes unreliable for temperatures of the order of the 
Fermi energy, which are typical values we are interested in; thus we have 
utilized the chemical potential provided by equation (|3.9|). 

Concerning the evaluation of the pair correlation function, the alternative 
method we have employed here (in analogy with the previous Section) con- 
sists in solving the Bethe-Goldstone equation for the relative motion of two 
interacting quarks. At T ^ the Pauli operator, which limits the available 
states in the formal development ( |2.1|) of the correlated wavefunction, will 
be modified by employing the smooth distribution function (|3.7|): 



0(k a - k F )6{k p - ki 



1-nUT) 



nr. 



XT) 



(3.11) 



As a result the effect of Pauli correlations is somewhat weakened and the 
influence of the medium can be expected to be, at a fixed density, less im- 
portant at finite temperature than at T = 0. 



9 



3.0 

<E> 

2.5 



2.0 



1.5 



1.0 



0.5 



0.0 

0.2 0.4 0.6 0.8 1 

P 

Figure 3: The ground state energy per particle at various (finite) tempera- 
tures, calculated in the independent pair approximation is shown 
as a function of the density: black dots correspond to T/Tc = 0.1, 
stars to T/T c = 0.5 and crosses to T/T c = 0.995. The dashed 
lines represent the corresponding free kinetic energy. 

Moreover, as already anticipated in the Introduction, as the temperature 
increases, one should also take into account some temperature dependence of 
the effective coupling, which we have assumed in the form suggested by ( |1.2|) . 
This potential has been employed in the evaluation of the two-body corre- 
lation function derived from the solution of the Bethe-Goldstone equation, 
where, as stated above, we also account for the temperature dependence of 
the occupation probability, according to ( |3.11| ). 

In addition, from the corresponding energy eigenvalues we have evalu- 
ated the ground state energy per particle as a function of the density, to be 
compared with the pure kinetic energy of the Fermi gas, and followed its 
evolution with increasing temperature (see Fig. 3). At high densities the sys- 
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tem is clearly a gas of non-interacting quarks; at low densities, instead, for 
small (and zero) temperatures the effect of pair correlations is clearly visible, 
but tends to fade away, as it is obvious from the temperature dependence 
of the quark-quark interaction, as the temperature approaches the critical 
value (here arbitrarily fixed "a priori" to coincide with twice the Fermi en- 
ergy at p = 0.5, T c = tt 2 /4). Notably the discontinuity in the derivative of 
the energy per particle, which is seen at T = and p ~ 0.2, is smoothed out 
at T ^ 0. 

Although, in the present treatment, we cannot identify a specific order 
parameter, which would allow to consider the phase transition from a ther- 
modynamical point of view, we have interpreted as a transition density, pc, 
the one where the medium induced correlations vanish: this quantity can 
be defined as the average value (with respect to the relative and total mo- 
mentum of a pair) of the difference between the matrix elements of the bare 
potential ( |1 . 1| ) and the G-matrix which one obtains from the solution of the 
Bethe-Goldstone equation: 



AU(p) = < k,K\V\k,K > - < k, K\G\k, K >. (3.12) 
In the above, according to the usual definition, 

l r R 

< k,K\G\k,K >= — / dxsin(kx)V(xU ok (x), (3.13) 
R Jo 

where both the unperturbed and the correlated wavefunctions are normalized 
over the above mentioned distance R and the integral over the center of mass 
coordinate is unity; the dependence upon K of the right hand side of ( 3.13 ) 



is implicit in the Bethe-Goldstone wavefunction through the action of the 
Pauli operator. 

The choice of AU as an "order parameter" is arbitrary, but it is closely 
related to the energy gap usually considered in the microscopic description of 
superconductivity; in that case a non-vanishing gap signals the existence of 
bound electron pairs which profoundly alter the global properties of the sys- 
tem. Here we assume ( |3.12j ) as a discriminant between a system of "hadrons" 
(bound pairs of quarks) and a weakly interacting Fermi gas of quarks. For 
each value of the temperature, the solution of the equation AU (p) = pro- 
vides the critical density for the deconfmement phase transition. 

The behaviour of pc with the temperature is shown in Fig. 4, which can 
be interpreted as the phase diagram for the model under investigation. For 
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Figure 4: The transition density versus temperature (in units of 7b), as it 
can be deduced by the quantity AU : the error bars are inherent 
to the averaging procedure employed in the determination of pc- 



a wide range of temperatures the critical density remains fairly constant, 
then it rapidly drops to zero as one approaches the critical temperature. It 
should be noticed the slight increase of pc in the low temperature regime: 
this outcome, somewhat opposite to the intuitive expectation, stems from 
the fact that a non-vanishing temperature weakens the Pauli blocking on 
the correlated wavefunction while the interaction strength is still practically 
unaffected. Thus for low temperatures the dominant effect is the appearance 
of the Fermi gas component at densities higher than for T = 0; as the 
temperature increases, however, this tendency is balanced and overcome by 
the temperature dependence of the interaction itself. 

Finally we have evaluated the pair correlation function, eq. (j2.7|) , by uti- 
lizing the finite temperature Bethe-Goldstone wavefunctions for the relative 
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Figure 5: The pair correlation function g(r) derived from the Bethe- 
Goldstone wave functions at T/Tc = 0.1 is displayed as a function 
of the relative interquark distance: the various curves are labelled 
by the value of the density. 



motion. A few examples of g(r) at various densities are illustrated in Figs. 5- 
7 for typical values of the temperature (in units of the above defined critical 
temperature). For low T (T = 0.1T C ) the same behaviour as for T = 
is observed, although, in agreement with the phase diagram of Fig. 4, the 
transition density appears to be somewhat larger than in the zero tempera- 
ture limit. As the temperature increases the binding effects become weaker 
and weaker and the correlation function rapidly approaches the free Fermi 
gas one, since the quark-quark model interaction becomes vanishingly small. 
Notice that the finite temperature washes out completely the quantum oscil- 
lations of the free correlation function. 
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Figure 6: The same as in Fig. 5, at T/T c = 0.5 

4 Remarks and conclusions 

In this work we have investigated some specific properties of a system of 
fermions, strongly interacting with a density dependent force, such as to 
provide bound pairs in the low density limit e a free Fermi gas at large 
densities. As already pointed out in a previous work fllf , this could bear 
some relevance in the investigation of the evolution of the primordial quark- 
gluon plasma toward the hadronic phase. However the evolution of the early 
Universe non only implies a decreasing local barionic density but also a rapid 
fall off of the temperature: the latter indeed should play a non-negligible role 
in the phase transition which lead to the formation of hadrons and nuclei. 

Here we have thus focussed the attention on the effects of temperature, 
starting from the hypothesis (supported by lattice QCD calculations) that 
with increasing temperature the strength of the confining potential becomes 
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Figure 7: The same as in Fig. 5, at T/T c = 0.95 



weaker and finally vanishes at some critical temperature. On this basis we 
have analyzed the interplay of density and temperature dependences of the 
quark-quark interaction: the pair correlation function shows indeed that, 
with increasing temperature, the transition density remains fairly stable until 
the critical temperature is approached and then rapidly drops to zero. 

The above findings show that the present approach is on the right path 
to achieve a microscopic description of the phase transition from hadrons to 
a quark-(gluon) plasma. However it can only be considered as a very prelim- 
inary stage: indeed the model for quark matter employed here suffers from 
many major limitations. Beside being one-dimensional and non-relativistic, 
it docs not embody internal degrees of freedom like spin, colour, etc. Also it 
does not contain explicit gluonic degrees of freedom, not even in the form of 
a perturbative one gluon exchange. 

Moreover, with respect to the present and future experiments of rela- 
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tivistic heavy ion collisions, one would also like to describe the dynamical 
evolution with time of the many-quark system, so that estimates of the 
probability for QGP formation and subsequent hadronization could be done. 
Future work is required in order to improve along these directions. 
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